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Using density functional theory (DFT), we analyze the influence of five classes of functional groups, as 
exemplified by NO2, OCH3, CH3, CCI3, and I, on the transport properties of a 1,4-benzenedithiolate (BDT) 
and 1 ,4-benzenediamine (BDA) molecular junction with gold electrodes. Our analysis demonstrates how ideas 
from functional group chemistry may be used to engineer a molecule's transport properties, as was shown 
experimentally and using a semiempirical model for BDA [Nano Lett. 7, 502 (2007)]. In particular, we show that 
the qualitative change in conductance due to a given functional group can be predicted from its known electronic 
effect (whether it is a/n donating/withdrawing). However, the influence of functional groups on a molecule's 
conductance is very weak, as was also found in the BDA experiments. The calculated DFT conductances for the 
BDA species are five times larger than the experimental values, but good agreement is obtained after correcting 
for self-interaction and image charge effects. 



One of the main advantages of molecular based electronics 
over present day semiconductor technology is the vast flex- 
ibility in design and functionality offered by the myriad of 
available molecules^. In order to fully exploit this flexibility 
for the synthesis of molecular devices, it is vital to establish 
simple guiding rules to estimate the effect that a given change 
in the molecule's composition or conformation has on its elec- 
trical properties. 

So far, most experimental and theoretical work on single- 
molecule conduction has focused on understanding the basic 
properties of individual junctions^ On the other hand, few 
studies hav e aim ed at describing general trends across molec- 
ular specieP^. 

It is well known that the chemistry of a functional parent or 
parent molecule may be altered in a consistent manner by the 
attachment of functional group^. Further, the influence of a 
particular functional group on a parent molecule's chemistry 
may be predicted qualitatively by considering the electronic 
effect of the functional group. 

We will demonstrate how such ideas may also be applied 
to describe the influence of a functional group on a parent 
molecule's transport properties. This was recently shown ex- 
perimentally for BDA, and the results were explained by scal- 



TABLE I: Functional group categorization by electronic effect with 
further examples ^ 





Functional 
Group 


Structure 


Electronic 
Effect 


Other 
Examples 


N02 


nitro 


A 


TT withdrawal 


CN, COR, 
SO3R 


0CH3 


methoxy 




TT donation 


NR2, OR 


CH3 


methyl 


a donation 


C2H5, C3H7 
alkyls 


CC13 


trichloromethyl 




cr withdrawal 


CF3, NR+ 


I 


iodo 


i 


cr withdrawal & 
TT donation 


F, CI, Br 
halogens 



ing the molecule's empirical conductance value by the cal- 
culated ionization potentials in the gas phase for the BDA 
species with two bonded Au atoms By employing these 
concepts from functional group chemistry, we may engineer 
the electronic transport properties of a known functional par- 
ent. This may be done using a functional group whose fron- 
tier orbitals have the same symmetry as the conducting or- 
bital of the parent molecule (a or tt). Groups with do- 
nating/withdrawing electronic effects may then be used to 
raise/lower the eigenenergy of the conducting orbital. 

In this paper, we analyze the influence of functional groups 
on the conductance of both a 1,4-benzenedithiolate (C6H4S2 
or BDT) and a 1,4-benzenediamine (C6H4(NH2)2 or BDA) 
molecular junction between gold contacts. We consider these 
molecules because BDT i s perhaps the best studied molecule 
for charge transporP^^, while the influence of functional 
groups on the BDA molecule has been recently studied experi- 
mentally In particular, we shall consider prototypical func- 
tional groups from each of the five electronic categories, as 
shown in Tablejl] These categories are: (1) withdrawal by con- 
jugation through the tt network (tt withdrawal), (2) donation 
by conjugation through the tt network (tt donation), (3) dona- 
tion by inductive effect through the a network (a donation), 
(4) withdrawal by inductive effect through the a network (a 
withdrawal), and (5) withdrawal by inductive effect through 
the a network with donation by conjugation through the tt net- 
work (a withdrawal & tt donation) ^ For each electronic effect 
categorization we have selected as respective examples a ni- 
tro group (NO2), a methoxy group (OCH3), a methyl group 
(CH3), a trichloromethyl group (CCI3), and an iodo group (I). 
Other examples of functional groups from each category are 
also provided in Table [l| 

All calculations have been performed using the SIESTA 
density functional theory (DFT) code^^ with the PBE ex- 
change correlation (xc) functionaP^, a double zeta polarized 
(DZP) basis set, and a mesh cutoff of 200 Ry. We first ad- 
sorbed the molecules on a five layer thick gold (111) slab and 
relaxed the molecule and the three outermost surface layers. 
The supercell contained 3x3 gold atoms in the surface plane, 
within which we used 2x2 /c-points. A pyramid of four gold 
atoms attached to the five layers of gold (111) was then intro- 
duced in the supercell to simulate the tip of a scanning tun- 
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FIG. 1 : (Color online). Schematics of (a) BDT species and (b) BDA species bonded between a gold (111) surface and tip, and their transmission 
functions versus energy E — Ep in eV, relative to the metal's Fermi energy Ep. 



neling microscope (STM). The pyramid was placed with its 
apex atom close to the unbound S/N atom of the molecule, 
and upon further relaxation the molecule bonds to the gold 
apex atom, as shown in Fig.[T] 

In this way we model a low temperature conductance mea- 
surement for a low molecular coverage on a gold (111) sur- 
face using an STM tip^^. Such a configuration may not accu- 
rately describe the results of break-junction experiments, es- 
pecially for the weakly bound BDA species However, as 
shown in Ref . JJ^ such a binding site yields a similar trans- 
mission function to more energetically favorable geometries. 
Thus, our chosen geometry should allow us to discern more 
broadly applicable trends in the electronic effects of the func- 
tional groups, which are our primary interest. 

We find that the BDT species (BDT, NO2BDT, OCH3BDT, 
CH3BDT, and I BDT) prefer the bridge site of the (111) sur- 
face and are rotated by approximately 30° to the surface nor- 
mal as shown in Fig.[T|a). On the other hand, we find that 
the BDA species (BDA, NO2BDA, OCH3BDA, CH3BDA, 
and I BDA) prefer the atop site and lie down at an angle of 
approximately 50° to the surface normal^^, as shown in Fig. 

Having chosen the contact geometries, we calculated 
the elastic transmission functions using the non-equilibrium 
Green's function (NEGF) formalism. The calculation proce- 
dure is equivalent to the one described in Ref. 17, except that 
the Hamiltonian was represented in the SIESTA atomic or- 
bital basis instead of Wannier functions. We note here that 
the DZP SIESTA basis set has recently been shown to yield 
transmission functions in quantitative agreement with plane 



TABLE II: Conductance G of BDT and BDA species between a gold 
(111) surface and tip. 



BDT Species 


G [2eV/i] 


BDA Species 


G [2eV/i] 


BDT 


0.28 


BDA 


0.024 


NO2BDT 


0.16 


NO2BDA 


0.024 


OCH3BDT 


0.32 


OCH3BDA 


0.026 


CH3BDT 


0.29 


CH3BDA 


0.024 


CCI3BDT 


0.25 


CCI3BDA 


0.022 


I BDT 


0.29 


I BDA 


0.023 



wave codes and maximally localized Wannier functionJ^. 
The transmission function was averaged over four /c-points 
in the surface plane, and three gold (111) layers were added 
on both sides of the molecule before the system was coupled 
to the bulk leads to ensure a smooth matching of the effective 
DFT potential. 

The calculated transmission functions are shown in Fig.[T] 
For both BDT and BDA the change in the transmission func- 
tion, T{E), near the Fermi level when the functional groups 
are attached, is strikingly small. 

The calculated conductances, G = GqT{Ef), where Go = 
2e^ /h, for the BDT and BDA species are shown in Table|ll| As 
is often the case for DFT transport calculations on molecular 
contacts, our results differ substantially from the experimental 
values of 0.011 Go obtained for BDT^\ and 0.0064 Go for 
BDA^^. We will address this issue in the last part of the paper. 
At this point, we observe that functional groups whose valence 
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FIG. 2: (Color online). Eigenenergies in eV and isosurfaces of ±0.02 
e/A^ of HOMOs for BDT species (squares) and BDA species (cir- 
cles), when adsorbed between a gold (111) surface and tip relative 
to the metal's Fermi level Ep (black dashed lines), and in the gas 
phase relative to vacuum (blue dotted lines). For BDA species, the 
experimental (diamond) and calculated (thin red solid line) gas phase 
ionization potential Iq with a SO shift of So ~ -2.7 eV, and the 
eigenenergies in the junction with a SO shift of So + ASo ~ -1.8 
eV (thick red solid line), are shown. 



orbitals are of tt symmetry seem to have the most influence on 
the molecule's conductance, with the tt withdrawing group 
lowering the conductance, and the tt donating group raising 
the conductance. 

For both BDT and BDA species, we found that transport 
primarily occurs through the HOMO, which in both cases is a 
TT bonding molecular orbitals, as shown in Fig.|2] In principle, 
the concept of molecular orbitals is somewhat artificial for a 
chemisorbed molecule, as hybridization effects will broaden 
the levels into resonances. One way of generalizing the con- 
cept is to consider i^moi, the projection of the Hamiltonian of 
the contacted system onto the subspace spanned by the ba- 
sis functions of the molecule. The eigenvectors of H^q\ can 
be regarded as molecular orbitals renormalized by the elec- 
trodes^^. For molecules which are not too strongly coupled to 
the electrodes, these renormalized molecular levels are easily 
identified with the levels of the free molecule. 

In Fig. [2] we compare the position of the HOMO level (rel- 
ative io Ef) of the contacted molecules with the HOMO level 
(relative to vacuum) of the free gas phase molecules. The 
variation in the HOMO position with the functional group 
correlates directly with the variation in the conductance (see 
Table [ll]) which shows that the current is indeed carried by 
the HOMO. One exception from this trend is NO2BDA for 
which the conductance is not lowered even though the HOMO 
is significantly down shifted. This is because the LUMO of 
NO2BDA lies so low that it also contributes to the charge 



transport. From Fig. [2] we can also see that each functional 
group produces a remarkably consistent shift of the conduct- 
ing orbital for BDA and BDT. This agrees with the basic 
premise of functional group chemistry, that a given functional 
group will alter different parent molecules in a consistent way. 
We also notice that the two tt groups (NO2 and OCH3) give the 
largest shifts of the eigenenergies, while the functional groups 
with cr symmetry have little influence. Moreover, the tt with- 
drawing functional group (NO2) delocalizes the HOMO and 
thereby lowers the eigenenergy, while the tt donating func- 
tional group (OCH3) confines the HOMO and thereby in- 
creases the eigenenergy. Thus the main ideas of functional 
group chemistry may be employed to qualitatively predict a 
functional group's influence on a given functional parent's 
conducting orbital(s). 

We stress that the qualitative effect of the functional groups 
on the HOMO of the contacted molecule, and thus its conduc- 
tance, roughly follows from the effect of the functional groups 
on the free molecule's HOMO, as was found experimentally 
for BDA^^. This is significant, as it suggests that the latter 
could be used as a simple descriptor to estimate the impact of 
a given functional group on the conductance, thereby allow- 
ing for an efficient screening of large numbers of functional 
groups. 

When we compare in Fig. [2] the magnitude of the shifts due 
to the functional groups for both BDT and BDA, we find sig- 
nificantly smaller effects for the contacted molecules than for 
the gas phase molecules. For example, the shift due to NO2 on 
BDT is approximately 0.4 eV in the gas phase, while it is less 
than 0.1 eV in the contact. This may be understood by rec- 
ognizing that the gold contacts act as electron sources/sinks 
for the parent molecule, counteracting the functional group's 
electronic influence. We thus believe that the weak influence 
of the functional groups on the conductance is a result of the 
self-consistent 'pinning' of the HOMO level which ensures 
the charge neutrality of the molecule. Such effects are indeed 
physical, but could be artificially enhanced by self-interaction 
errors in the PBE xc-functional: For partially occupied molec- 
ular orbitals, an incomplete cancellation of the Coulomb self 
interaction by the xc-functional will artificially raise (lower) 
the energy of that orbital when charge is added (removed) and 
this will enhance the 'level pinning'. 

Self interaction errors also contribute to the well-known un- 
derestimation of band gaps by DFT calculations^^. Recen tly, 
an atomic self interaction correction (ASIC) schemJ^^^ has 
been proposed as a simple cure to this problem for molecular 
contacts. In general, however, image charges formed in the 
metallic electrodes when electrons are added to or removed 
from the molecule also renormalize the molecular levelP, 
and this effect is not captured by the SIC. 



In Table [ni| we compare our DFT calculated conductances 
for the BDA species with corresponding experimental val- 
ues The calculated numbers are roughly five times larger 
than the experimental ones due to the above mentioned de- 
ficiencies of the DFT approach. It is, however, interesting to 
notice that: (i) the weak effect of functional groups is found in 
both data sets (ii) the qualitative effect of the functional groups 
on the conductance (increase/decrease) is reproduced by cal- 
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TABLE III: (Color online). Calculated conductance Ccaic, SO con- 
ductance Gso^^ with shift Eo — AEo ~-1.8 eV, and experimental 
conductance Cexp^*^ for BDA species between a gold (111) surface 
and tip. 



BDA Species 


Gcaic [2eV/i] 


Gso [2eV/i] 


Gexp [2eV/i] 


BDA 


0.024 


0.0035 


0.0064^ 


NO2BDA 


0.024 


0.0032 




OCH3BDA 


0.026 


0.0035 


0.0069^ 


CH3BDA 


0.024 


0.0034 


0.0064^ 


CCI3BDA 


0.022 


0.0034 




CI BDA 


0.015 


0.0025 


0.0060^ 


I BDA 


0.023 


0.0036 




^Ref.fTol 



culations. In the second column of Table III we present the 
results of our conductance calculations after a 'scissors opera- 
tor' (SO), which corrects for self-interaction errors and image 
charge effects, has been applied to the molecule's spectrum. 

To construct the SO we follow a recent work by Quek et 
ah who applied the scheme to a BDA-gold junction In this 
method, the underestimation of the DFT HOMO and LUMO 
are corrected by shifting the occupied/unoccupied orbitals of 
the contacted molecule byi;Q'^ = ±i;o = T(^o + ^homo), 
where Shomo is the DFT HOMO level for the free molecule 
and /o the ionization potential^^ . We also calculate the later 
in DFT from /q = ~ ^q=o. where is the to- 

tal energy of the molecule with charge +e in the gas phase, 
while Eq^o is the total energy in the gas phase of the neu- 
tral molecule. For BDA without functional groups we find 
/o ^ 6.73 eV, which agrees quantitatively with the experi- 
mental value of 6.83 eV^. We found only small variations in 
So for the different BDA species and have used Sq ^ -2.7 eV 
for all molecules, as indicated in Fig. [2] 

To estimate the shift of the molecular levels by image 
charge effects, we first calculate the charge distribution af- 
ter removing an electron from the free molecule. A Mul- 
liken analysis is then used to approximate this continuous 
charge distribution by point charges located at the atoms of the 
molecule. We model the gold contacts as two perfectly con- 
ducting surfaces separated by 13.2 A, which corresponds to 
the vacuum separation between the opposing gold (111) sur- 
faces. We then obtain an image potential AHq ^ -0-9 eV for 
all of the BDA species considered^"^. For the BDA molecular 
junctions, we thus shift all occupied/unoccupied orbitals (ob- 
tained by diagonalizing i^moi) by Sq'^ — ASq'^ = ^(So — 
ASq) ^ eV as indicated in Fig. [2] Calculating the 

conductance using the renormalized Hamiltonian yields the 
values shown in Table HHl We find that this 'ad hoc' correc- 
tion produces results approximately half those obtained from 
break-junction experiments on the BDA species^^. This is 
most probably due to our choice of contact geometry, as our 
calculated conductance for BDA of 0.024 is about half the 
average DFT value calculated over 15 different break-junc- 
tion contact geometries of 0.046 Go We stress, how- 
ever, that the use of this SO can only be justified for weakly 



coupled molecules with HOMO/LUMO levels well separated 
from Ep. In Fig.[2]we see this is the case for the BDA species, 
since the HOMOs are localized near the amine contact groups, 
providing a poor overlap with the gold contact orbitals. How- 
ever, for the BDT species, whose HOMOs are rather diffuse 
around the sulfur atoms, there is a strong overlap with the gold 
contact orbitals, so that application of the SO is not justified. 

In conclusion, we found that the ideas of functional group 
chemistry may be applied to qualitatively predict the influence 
of a functional group on the electronic structure of a parent 
molecule, as has been found for BDA both experimentally and 
using a semiempirical model However, we also found that 
functional groups have a very weak influence on a molecule's 
conductance in agreement with recent experiments^^. The rea- 
son for the weak influence is that charge neutrality pins the 
HOMO/LUMO molecular levels, making it difficult to shift 
them relative io Ep. By applying multiple functional groups 
to the same parent molecule it may be possible to obtain a 
stronger influence, as was found for BDA^^. By employing 
a scissors operator correcting for self-interaction errors and 
image charge effects, we obtained qualitative agreement with 
experimental conductance values for the BDA species. Our 
results suggest that effective 'switching' of a molecule's con- 
ductance may require a more direct change in the strength or 
geometry of the molecule's contacts, to overcome the 'level 
pinning' of the metal. 

We thank J. K. N0rskov, S. Dobrin, and M. Strange for 
useful discussions. The authors acknowledge support from 
the Danish Center for Scientific Computing through grant No. 
HDW-1 103-06. The Center for Atomic-scale Materials De- 
sign (CAMD) is sponsored by the Lundbeck Foundation. 



APPENDIX A: SCISSORS OPERATOR (SO) 

The SO for a free molecule in the gas phase So, is given by 

So = -(s^HOMO + ^0), (Al) 

where ^homo is the DFT HOMO level of the free molecule 
in the gas phase, and /o is its first ionization potential, which 
is defined as 



/o = ^< 



E^ 



q=0' 



(A2) 



Here £^g=+i is the total gas phase energy calculated for a 
molecule with charge q = +e, while Eq=o is the total energy 
in the gas phase for the neutral molecule. 

The occupied/unoccupied orbitals of the contacted 
molecule are shifted by the same amount Hq'^ = iSo for 
each orbital. This is reasonable for orbitals other than the 
HOMO and LUMO, since these orbitals are sufficiently 
far from the Fermi level Ep to have little influence on the 
contacted molecule's conductance. 

However, this is not the case for the LUMO orbital, espe- 
cially for NO2 BDA. In this case, a fully rigorous calculation 
would shift the DFT LUMO level by -(slumo +^ea), where 
s^LUMO is the DFT LUMO level for the free molecule and E^a 
the electron affinity. The later may be calculated similarly to 



the ionization potential in DFT from Eea = E^ 



q=0 



Ea 
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where Eq^_i is the total energy of the molecule with charge 
— e in the gas phase. However, as discussed in Ref. [TSl 
when such a shift was calculated for the BDA LUMO, it was 
found to be approximately equal and opposite to that for the 



HOMO. For these reasons, we use Sn 



±So for all the 



occupied/unoccupied orbitals. 

For a molecule adsorbed near a metal surface, the SO is 
So —A So, where A So is the charged molecule's image poten- 
tial arising from the metal surface. To estimate the charge dis- 
tribution on the molecule p(r), we employ the gas phase Mul- 
liken analysis from a charged spin polarized calculation, and 



use the distribution of unpaired spin up charges. Table IV a) 



and IV ^b) show the distribution of charges on each atom in 
(a) 1,4-benzenediamine (BDA) and (b) 3-nitro-l,4-benzene- 
diamine (NO2BDA) respectively. 

The total potential energy gain by moving the charged 
molecule to the metal surface is 

ASo = /d^rp(r)$ind(r) + ^ /d^rpind(r)^ind(r), (A3) 



where the first term is the self energy of the molecule's charge 
distribution p(r) due to the electrostatic potential ^ind(r) in- 
duced on the metal surface. 



*^>ind(r) 



and the second term is the Hartree energy for the charge 
distribution Pind(r) induced on the metal surface. Approx- 
imating our gold surface by a perfect conductor, we find 



TABLE IV: Mulliken analysis of charge distribution on (a) BDA and 
(b) NO 2 BDA in the gas phase with total charge q = +e. 



(a) 1,4-benzenediamine 



Qi [e] Xi [A] Vi [A] Zi [A] 




+0.233 


2.508 


2.532 


4.225 


+0.153 


3.435 


2.297 


5.257 


+0.057 


3.770 


0.972 


5.657 


+0.057 


4.027 


3.390 


5.954 


+0.057 


4.766 


0.754 


6.635 


+0.057 


5.028 


3.165 


6.916 


+0.153 


5.438 


1.846 


7.243 


+0.233 


6.553 


1.663 


8.141 



(b) 3 -nitro, 1,4-benzenediamine Qi [e] Xi [A] yi [A] Zi [A] 



+0.237 


2.591 


2.540 


4.187 


+0.142 


3.504 


2.302 


5.220 


+0.007 


4.089 


3.397 


5.932 


+0.100 


3.834 


0.989 


5.626 


+0.076 


5.071 


3.177 


6.909 


+0.070 


4.822 


0.773 


6.626 


+0.015 


4.361 


-1.509 


6.479 


+0.082 


5.506 


1.865 


7.258 


+0.263 


6.601 


1.694 


8.149 


+0.008 


5.979 


-0.878 


7.847 



the induced charge distribution is restricted to the surface, 
Pind(r) = cr(r^)(5(x), where cr(r^) is the induced surface 
charge distribution, and = yey + ze^ are the components 
of r in the ?/z-plane, so that 



.(A5) 



ASo = jd\p(v)^,^^(v) + iy"d2r^a(r^)$ind(r^) 

Using the boundary condition that the total electrostatic po- 
tential is zero on the metal surface, we find ^ind(r^)|ic=o = 
— ^moi(r^) |x=o, where <^>moi(r) is the Coulomb potential due 
to the charged molecule. We then find 

AEo = y"dV(r)*ind(r)- d\^a{v^)^rno\{v^) ,(A6) 
= y"d^rp(r)<E>ind(r) - ^y"d3rpind(r)$moi(r), (A7) 

(A8) 



Since we are modeling the molecule's charge distribution 
as a series of N point charges at locations = {xi^yi^Zi} 
with charges Qi, so that 



N 



(A4) p(r) = ^Q,(5(r-ri), 



(A9) 



and we will only be evaluating the induced potential outside 
the metal, we may employ the method of image charges, so 
that 



(AlO) 



where we have replaced the in duced charge density by the im- 
age charge density Pind(i') in (AlO). We may then write 



N 

Pimg(r) = ^Qjp[4W' 
so that 

AT 



(All) 



(A12) 



- + 



- + 



- + 



-2L-X -2L+X 



2L-X 2L+X 



4L-x4L+x 



FIG. 3: Locations of image charges for two perfectly conducting 
surfaces at x = and x — L due to a positive charge at +x, where 
L — 13.2 A for our system. 



6 



where Pimg(^) image charge density in the metal due to 
a point charge of 1 e located at vj . 

If the metal our molecule is adsorbed upon consists of a 
single perfect conducting surface located at x = 0, we may 
model the charge distribution induced in the metal as a series 
of image charges, so that 



-S{x + Xj)S{r^ 



1 ^ 



QiQj 



{xi + xj)^ + Rl 



(A13) 



(A14) 



where Ra = llr^ 



However, if we now introduce a second surface located at 
X = L, an infinite series of image charges is produced, anal- 
ogous to a series of repeated reflections between two mirrors. 
We then find 



-S{t^ - rf) ^ [S{x + Xj + 2(n - 1)L) + S{x + Xj - 2nL) + 5{x - Xj + 2nL) -^5{x + Xj - 2nL)] , (A15) 



n=l 
1 



{xi + x^- - 2nL)2 + /(^ 



•2(n-l)L)2+M. 



^{xi-Xj^2nLY^R}. ^{xi - Xj + 2nLf + 



(A16) 



as shown in Fig. [3] We shall take the separation between the 
opposing gold (111) surfaces of 13.2 A as our L value. It 
should be noted that the image charges produced in such a 



manner are not physical charges, but are employed as mathe- 
matical tools to ensure that the total electrostatic potential is 
zero on both perfect conducting surfaces. 
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